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ABSTRACT 

The state after virialization of a small-to-intermediate A^-body system depends on its ini¬ 
tial conditions; in particular, systems that are, initially, dynamically “cool” (virial ratios 
Q = 2T/|n| below ^ 0.3) relax violently in few crossing times. This leads to a metastable 
system (virial ratio ~ 1) which carries a clear signature of mass segregation much before the 
gentle 2—body relaxation time scale. This result is obtained by a set of high precision N- 
body simulations of isolated clusters composed of stars of two different masses (in the ratio 
mh/mi = 2), and is confirmed also in presence of a massive central object (simulating a black 
hole of stellar size). We point out that this (quick) mass segregation occurs in two phases: the 
first one shows up in clumps originated by sub-fragmentation before the deep overall collapse; 
this segregation is erased during the deep collapse to re-emerge, abruptly, during the second 
phase that occurs after the first bounce of the system. This way to segregate masses, actual 
result of a violent relaxation, is an interesting feature also on the astronomical-observational 
side. In those stellar systems that start their dynamical evolution from cool conditions, this 
kind of mass segregation adds to the sequent, slow, secular segregation as induced by 2- and 
3- body encounters. 
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1 INTRODUCTION 

Astronomical observations indicate that several stellar systems 


(from young and very young open star 
of galaxies) show a certain degree of se 
minous and massive components in their 

clusters to rich clusters 
gregation of the most lu- 

inner regions i Bontemps 

^rari([20T0ll,|Kirk et al.| j20T4),|Er et al. 

^2013|l,|Gouliermis et al. 

(|2006k IRaboud & MermilliodU19981l, | 

laboudl (|1999k |Littlefair 


|et al.| j20^). For instance, the Orion Nebula Cluster (ONC) is 
mass segregated down to about 5 Mq despite its young age, es¬ 


timated to be less than 2 Myr ( [Hillenbrand & Hart mann| ( |1998| l, 
I Allison et al.lpOO^ ). The dominant cause of the rapid mass segre¬ 
gation process, for such systems, is still under dehate. A dynamical 
origin is commonly excluded; in fact, many of the systems which 
show a mass segregation are much younger than their two-body re¬ 
laxation time, which is usually considered as the time needed to 
segregate masses. Specifically, the relaxation time is defined as 



where v is the typical velocity of a star in the system, and 
D [(Auy)] is one of the three, independent, diffusion coefficients 
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which derive from the Fokker-Planck equation. For simplicity, 
if we assume that the velocity distribution of the field stars is 
Maxwellian, with dispersion a, it is possible to obtain an explicit 
expression for D [(Auy)] and Eq. 1 becomes 


irel 
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where p is the mean mass density of the field (target) stars, 
m is the mean mass. In A is the usual Coulomb logarithm, X = 
w/(\/2cr) and G{X) is the function 


G(^) = ^ 


2 v2 

erf(X) - —Xe~^ 

■s/Tl 


(3) 


Equation is valid whether the initial conditions are, in¬ 
deed, not too far from virial equilibrium (see for example |Binney| 
|& Tremaine|p008) l). This is not always a correct assumption, es¬ 
pecially if we refer to the study of the early dynamical evolution 
of stellar systems whose stars form in sub-structured, clumpy re¬ 
gions in sub-virial conditions (|De March! et al.| j2013f,|Schmeja &| 
|KlessM 1 200^ Adams et al.|j2006|). It haTbeen already shown 


|Farouki & Salpeter 1 1982 1 , Allison et al.H2009') , [McMillan et al.| 
2007^ ), through A-body simulations, that, if the initial state of 
a stellar system is out of equilibrium and the initial spatial dis¬ 
tribution of its stars is not homogeneous, a significant degree of 
mass segregation can emerge on very short time-scales, signifi- 
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cantly shorter than the two-body relaxation time. | Allison et al.| 
( |2009^ argue that the main mechanism that leads to mass segrega¬ 
tion on a short time is related to the short interval of time in which 
the violent collapse creates a very dense core; typically, it contains 
about half the mass of the stellar system in a radius of about one 
tenth of its initial size. They show that the time to segregate masses 
down to 4 — 5 Mq, in the dense core, is comparable to its liv¬ 
ing time (approximatively 0.1 Myr). This is enough to justify also 
the degree of mass segregation observed in some astrophysical sys¬ 
tems, such as the cited ONC ( [Allison et al.|j200^ ). On the other 
hand, [Bonnell & Davie^ ( |1998| ( excluded that the observed mass 
segregation in young stellar clusters could be due to a violent dy¬ 
namical evolution, introducing, rather, the hypothesis of an in situ 
formation of the most massive and luminous stars. In their work, 
they investigated the dynamical evolution of both spherical stel¬ 
lar systems, initially in virial equilibrium, and non-spherical stellar 
systems in sub-virial conditions. They found that the time-scale for 
mass segregation was largely unaffected by differences in the initial 
phase-space distribution of the stars. 

A recent paper by jCaputo et al.j j2014| l investigated the role 
of initial sub-virial conditions on the final state of a set of stellar 
system on a range of N from 2,048 to 131, 072 and intermediate 
mass black hole therein and found, by a matching of results with 
some observational data, that the initial values of the virial ratio 
Q = 2T /|fl| should be in the range 0.36-0.50. They also found an 
enhanced mass segregation in initially cooler systems but they give 
no evidence of sub-clustering before the collapse. 

In the light of the described scenario, it is clear that it is dif¬ 
ficult to discriminate between the two formulated interpretations to 
explain the phenomenon of the rapid mass segregation in young 
star clusters. Moreover, an additional problem comes from the ob¬ 
servational side, because very young star clusters, still sites of star 
formation processes, are very arduous to observe because they are 
embedded in high density gas clouds. Therefore, it is hard to prove 
if the most massive stars form, indeed, very close to the innermost 
regions or, on the contrary, they form on a large spatial scale and 
segregate later. 

Another dynamical mechanism, not deeply investigated yet 
(even if already highlighted by [Aarseth et al.j jl 988^ , jMcMilla n] 
jet al.| ( [2012| () and that might play an important role in segregating 
masses, is the initial rapid fragmentation of a stellar systems whose 
stars are initially distributed homogeneously with very small ini¬ 
tial velocities. [McMillan et al.| ( |2007| l and [McMillan et ^j2012^ 
stressed that several sub-systems form during the collapse phase; 
in particular, they show, just before the bounce of the system, a 
certain degree of mass segregation which is preserved after the se¬ 
quent collapse and rebound. It is important to stress that this kind 
of segregation occurs before the formation of the short-living core 
described in [Allison et al.[ ( [2009[ l. 

In this context, we are forced to conclude that the theoreti¬ 
cal scenario is still unclear and deserves investigation. Moreover, 
the real situation is made even more complicated by the presence 
of a certain number of primordial binaries and/or of black holes, 
and/or of a background gas (not considered in A^-body simulations, 
so far). In particular, in this paper we study the effects of the violent 
collapse of an A-body self-gravitating system starting from initial 
conditions corresponding to a virial ratio 0 ^ Q 5/ 1 with stars uni¬ 
formly distributed in space. We also check the role played by the 
presence of both a residual gas after star formation and of a stel¬ 
lar mass black hole on mass segregation and on the stars’ velocity 
dispersion. 

In all our A-body simulations, privileging the clarity of re¬ 


sults respect to their generality, we consider the simple case of a 
bimodal mass spectrum with bodies initially distributed randomly 
in a sphere of unitary initial radius R. 

The organization of the paper is: in Sect. we describe the 
models of stellar systems we adopted and give a description of both 
the software and hardware resources used to follow their dynami¬ 
cal evolution; in Sect. we present and discuss the results with 
attention to both the physics of violently relaxing, intermediate A, 
systems and to the possible comparison with observational data of 
real clusters. Finally, in Sect. we list our conclusions and outline 
the developments needed to get a better insight into the evolution 
of open clusters emerging from their mother proto-clouds. 


2 MODEL 

2.1 Stellar systems 


We performed a large set of direct A-body simulations of young 
star clusters with 128 ^ A ^ 1,024. We adopted a bimodal mass 
spectrum: “light” stars (each with mass mi and total number Ni) 
and “heavy” stars (each with mass mn and total number Nh) such 
that mh = 2mi and Nh = A;. The initial spatial coordinates of 
bodies have been sampled from a uniform distribution in a sphere of 
radius i? = 1. The initial velocities have random orientations and, 
in order to determine their absolute values, we assume a certain 
value of the initial virial ratio Q, defined as 


Q = 


2T 


(4) 


where T is the kinetic energy and the potential energy of the A- 
body system. We vary Q from Q = 0, corresponding to the most 
violent (cold) collapse, up to Q = 1, which means virial equilib¬ 
rium, at steps of 0.1. Small values of Q (“cool” systems) seem to 
be suggested by observations of young stellar systems where stars 
appear to be clustered in sub-virial conditions (see, for example, 
[Adams et al.[ ( [2006t and [Hillenbrand & Hartmann| ( [1998[ l). We also 
investigated the gravitational role of a residual gas, modelled as 
an analytical contribution to the accelerations of stars. This field 
is represented, at any time t, by a time varying Plummer potential 
(|K ummer HMD) 


<E>P (r;/) 


GMg 

Vr^+^c (i) 


(5) 


where (t) is the core radius, assumed time-dependent, and Mg is 
the total (constant) gas mass. Another ingredient of our simulations 
is the inclusion of a particle with a mass significantly larger than the 
others. This object may be considered as a stellar black hole, whose 
mass is indicated as msH- The black hole mass was assumed to be 
mBH = 25mh = 50mi. The total mass of the star cluster system. 
Mg, is assumed as unit of mass, that is 


Mg = Nhmh -\- Nimi + Mg + msH = M* -b Mg = 1 (6) 


where M* refers to the total mass in stars. The mass of the gas the 
stellar system is embedded in is given by assuming a value for the 
star formation efficiency parameter, that is 


M, _ M, 
Mg Mg Mi, 


(7) 


Once given 5, the stellar mass of the cluster is M* = SMg and 
Mg = Mg{l — S). We considered two values of <S, 5 = 0.3 and 
<S = 1 (no gas). Because we are interested in the emerging state 
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of young open clusters, we evolved our stellar systems up to a rel¬ 
atively short time (< 5Myr), therefore we neglected, in a first ap¬ 
proximation, the effects of stellar evolution. For each set of values 
of the free parameters we performed 30 runs (> 5,000 simulations, 
in total) to give a statistical significance to our results. 

The results of our simulations will he presented using the 
crossing time of the system as time unit, that means for the gravi¬ 
tational constant the value G = 1. The velocity units are obtained 
consequently. 


2.2 The A'^-body code 


For the purposes of this paper, we used our highly parallel Al-body 
code HiGPUs ( [Capuzzo-Dolcetta et al.| ( |^13^ ) running on our pri¬ 
vate machine containing a Central Processing Unit (CPU) Intel i7 
950 and 2 nVIDIA Tesla C2050 (Fermi) Graphics Processing Units 
(GPUs). Although the GPUs used are not the best in terms of cost 
and computing capability, being old generation cards, thanks to the 
high single core working frequency they perform well in regime 
of weak load, i.e. using a number of particles < 1,024 jCapuzzo-| 
|Dolcetta & Spera|2013^. HiGPUs is base d on a Hermite 6 th order 
integration scheme i Nitadori & Makino||2008| > parallelized using 
CUDA (or OpenCL) plus OpenMP and MPI to guarantee maxi¬ 
mum performance when a hybrid (CPUs-l-GPUs) computer is used. 
The integration algorithm is implemented using the technique of 
block time steps ( |Aarseth|2003| > where the particle time steps are 
determined using the following formula: 


Afi = 


at + Q 2 


^(i)\ /^(i) 


A' 


1/3- 


( 8 ) 


where 

^(s) _ /I (s —1)1 I (s + 1) I I I (s)l^ /pvN 

A-l=y\a.\ l + r* I 

and is the ^-th time derivative of the acceleration of the i-th 
particle. Equation [^represents a weighted mean (with coefficients 
ai and 02 ) between the Aarseth criterion for the 4th order Hermite 
integrator (with accuracy parameter 774 ) jAarseth| ( |2003^ ) and the 
generalized Aarseth criterion for the 6 th order scheme (with accu¬ 
racy parameter tyg) ( |Nitadori & Makino]j2008| >). The combination 
with Q?! = a 2 = I has been found to be more stable, for the 6 th 
order method, than the two criteria used independently, providing a 
better total energy conservation and avoiding time steps either too 
large or too small. 


3 RESULTS 

Because no significant differences were found when changing the 
number of stars in the range studied (128 ^ A 1, 024), in the 
following we will discuss the results only for N = 1, 024 that we 
choose as reference case. One good indicator to quantify the level 
of mass segregation is the ratio of the lagrangian radii of heavy par¬ 
ticles to those of light particles, so that values significantly greater 
than 1 indicate the presence of mass segregation. Of course, the use 
of lagrangian radii ratios as mass segregation indicators is reliable 
when the system has a spherical symmetry. Nevertheless, before 
the first bounce of the system, we observe a fragmentation of the 
system into various clumps which are not spherical. In such case, a 
more reliable evaluation of the presence of mass segregation comes 
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Time 


Figure 1. The ratio between and A; for several values of the initial 
virial ratio Q, in Case 1. Values of Q between 0.5 and 0.9 have also been 
considered, but they are not shown here for the sake of a clear representation 
of the results. Each curve represents an average value of Af^/Ai and the 
error (standard deviation) is represented by the semi-transparent area. 


from the so called minimum spanning tree method (MST|^ devel¬ 
oped by |Allison et al.| j2009| l. Given a sub-set Km of m points 
belonging to a system composed of m' > m bodies, the degree of 
mass segregation established for that sample, Ax„, is defined as 

A {tnorm'} . tTnorm /1 r»\ 

Km = - ± - ( 16 ) 

iKm l-Km 

where is the MST for the sample Km, (Inarm) is the aver¬ 
age MST for m randomly selected stars in the whole system and 
(Tnorm is its associated standard deviation. In each run, we calcu¬ 
late (Inarm) usiug 200 different sub-sets of points. From the defi¬ 
nition given in Eq.[^it follows that the generic sample K is mass 
segregated if Ak^ is significantly greater than 1 . 

When we use the MST method to investigate the distribution 
of masses in our simulations, we need to remove escapers. This re¬ 
moval process is important because, after the first bounce, a signifi¬ 
cant amount of mass (about 20 % of the total mass) is lost; therefore, 
a single, remote, star may alter significantly the length of the span¬ 
ning tree of a specific population. To identify escapers correctly, the 
best criteria are those based upon energy (which should be positive) 
and distance from the densest region of the system (the core). 


3.1 Case 1: systems without gas and without a central black 
hole 

The simplest case we studied concerns systems with neither gas nor 
a central black hole. To give statistical significance to our results, 
we generated a set of 30 different initial configurations to evolve, by 
simply changing the seed of the Mersenne Twister random number 
generator jMatsumoto & Nishimura| ( |1998^ ) that provides initial 
conditions in the phase space. Here we show the results obtained 
from these simulations. 

Fig.[T]shows the ratio between the averaged value of A for the 
heavy (Ah) and for the light (A;) stars as a function of time. The 
first important thing to stress is that the degree of mass segregation 


^ The minimum spanning tree is the shortest path length which connects a 
certain number of points without forming closed loops. 
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Figure 2. A detail of Fig.[T] zooming in the interval of time between 0 and 
2.5. 

depends strictly upon the initial state of the system: the farther from 
equilibrium, the higher and the quicker the degree of the resulting 
mass segregation on both short and long time scales is. Another 
important result is that we obtain a significant degree of mass seg¬ 
regation starting from homogeneous and smooth initial conditions. 
For the cases of most violent collapses (Q = 0 and Q = 0.1) the 
system gets to a saturation of the degree of mass segregation around 
12 units of time, while, for larger values of the initial virial ratio, 
the process of mass segregation continues at an approximatively 
constant rate. 

Fig.|^shows a detail of the Figj^ in particular, it refers to the 
initial evolution of the system from f = 0 to just after the bounce. 
This figure makes clear that there is a rapid increase of mass seg¬ 
regation up to f ~ 1; then, an inverse trend is observed for a brief 
time and, finally, mass segregation starts again with about the same 
rate and efficiency, as it was before. The series of snapshots in Fig. 
[^give a visual sketch of what explained above. 

The rapid increase of mass segregation before the bounce is 
in perfect agreement with the mass segregation observed in the 
sub-clumps that form during the first system collapse. Neverthe¬ 
less, when the merging process occurs, the mass segregation is not 
completely preserved, and Fig|^shows an inverse trend, indeed. Af¬ 
ter the merging process, mass segregation continues with the same 
efficiency as before the bounce but, this time, no sub-structures ap¬ 
pear; therefore, mass segregation continues inside the dense core 
formed just after the bounce. It is thus clear that sub-clumps are not 
the unique cause of the observed mass segregation on short time 
scales, and the same can be said for the short-living but very dense 
core. Actually, both the two phenomena play a role because the 
first one acts on a very short time-scales, f < 1, while the second 
one is responsible of the long lived mass segregation for f > is. 
Times around the collapse phase (f ~ fs) correspond to a transi¬ 
tion between the two regimes. Moreover, this situation is evident for 
“cool” systems (Q < 0.3); for Q > 0.3, sub-clumps do not form at 
all, although a degree of mass segregation significantly greater than 
that in the equilibrium case (Q = 1) is established. In these “warm” 
cases (Q > 0.3), the only mechanism active in the rapid and sec¬ 
ular mass segregation is the dynamical evolution of the dense core 
that forms after the collapse. 

For completeness, we report in Fig. I^the results obtained us¬ 
ing another indicator of mass segregation, that is the ratio between 
some lagrangian radii of light and heavy stars. As we said, the ratio 



Figure 3. Projection on one coordinate plane (x-y) of one of the simulated 
clusters with initial virial ratio Q = 0 at four different times. Panel (a) 
shows the homogeneous initial distribution, panel (b) the formation of sev¬ 
eral sub-clumps (t ~ 1), the most evident of which are circled, panel (c) 
refers to the state of maximum compression of the stellar system (t ~ ts), 
while panel (d) is a view of the cluster after the bounce (t ~ 2). 
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Figure 4. Time evolution of the ratio of the lagrangian radius (that contain¬ 
ing 50% of the mass) of light to that of heavy stars in the evolving system, 
in Case 1. Different curves are for different values of the initial virial ratio 
Q of the system. The pattern area represents the error associated with each 
curve. 


of lagrangian radii is as better as mass segregation indicator as the 
studied system maintain spherical symmetry. In our case, this is true 
at times greater than ~ is. In other words, we use lagrangian radii 
to point out the presence of a long lived mass segregation whose 
efficiency depends on the above described phenomena which occur 
at f < fs. 

The phenomenon of mass segregation is particularly evident, 
in Fig. 1^ for (5 = 0. On average, in this case, the typical size 
of the spatial region occupied by light stars becomes, in 6 crossing 
times, almost twice larger than the heavy particles one. On the other 
side, the case Q = 1 shows an almost flat behaviour. This confirms 
again, the strong dependence of the efficiency of mass segregation 
upon the initial dynamical condition of the system. 


3.1.1 Some analytical considerations 


It is known ( |Spitzer|I969| > that the segregation time of a sub pop¬ 
ulation of stars of mass M in a system where (m) is the average 
mass is 


(tTL/ 

tseg{M) ~ ^trel. (H) 

Nevertheless, if the system is far from equilibrium, as it is in our 
cases, further simplifications of the formula for the relaxation time 
(see Eq.|^ are not allowed. 

An expression which does not assume virial equilibrium, is 
the following 


trel — 0-^4 „ / \ 1 A 

{m) p log A 

while assuming virial equilibrium we get 


( 12 ) 


^ N R 
8 log N a 


(13) 


In our simulations, with = 1, 024, the dense core contains 
N ~ 350 stars and has a size i? ~ 0.15 pc. The velocity dispersion 
is (7 ~ 3 km/s and the average mass (m) ~ 1 M©. These values 
give trel = 0.3 Myr when using Eq. [T^ which means that the 
system evolves fast and can segregate masses, down to 1.5 M©, on 


Figure 5. The spanning trees ratio as a function of time, in Case 2. Values of 
Q label the curves. Values of Q between 0.5 and 1 have been studied, too, 
but displayed here for the sake of more clear representation of the results. 
Each curve represents an average value of A^/Ai, taking into account the 
results obtained from the single runs, while the error (standard deviation) is 
represented by the semi-transparent area. 


a very short time scale (the very dense core lives, in our simulations, 
for about 0.2 Myr). On the other hand, using Eq.[T^ we get a value 
of about 3 Myr for the relaxation time, therefore we have a system 
that can segregate masses down to 15 Mq only. This is just one 
proof of the inconsistency and inapplicability of both formulas |13| 
and|12|to a situation which is very far from virial equilibrium. 


3.2 Case 2: systems containing a central black hole without 
gas 

As we can see in Fig.|^ the presence of a central massive particle, 
whose mass is 25 times that of a generic “heavy” star, tends to 
compact the curves of the spanning trees ratios with respect to the 
case shown in Fig.[^ Actually, the process of mass segregation is 
less efficient but it is still evident and strongly dependent on the 
violence of the collapse. 


3.3 Case 3: systems containing gas without a central black 
hole 

The majority of young and very young star clusters are still em¬ 
bedded in their proto-cloud. This residual gas, when sufficiently 
abundant, acts as a background potential which can affect the dy¬ 
namical evolution of the stellar cluster. In this work we modelled 
the presence of a gas as an analytical potential to add to the pair 
gravitational interaction of the bodies in the system. This does not 
represent very correctly the real astrophysical situation but, using 
our simple model, we can at least account at order of magnitude 
accuracy the gaseous phase contribution. The most evident effect 
of the addition of the analytic gravitational potential is smoothing 
out the 2-body encounters, decreasing the efficiency of mass segre¬ 
gation on both small and longer time scales. 

Figure is the analogous of Fig. [T] with the inclusion of a 
stationary gaseous background. As we can see, the results for the 
case of (5 = 0 are very similar to those obtained for the same case 
but without the presence of gas; it is still evident the formation of 
sub-clumps and a quick mass segregation, before the bounce, as 
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Figure 6. The spanning trees ratio as a function of time, in a case analogous 
to Case 1 but with a relic, static gas. Values of Q between 0.4 and 0.9, 
both taken into account in our simulations, are not shown here to have a 
clearer representation of the results. Each curve represents an average value 
of A/f /A/,, taking into account the results obtained from the single runs, 
while the error (standard deviation) is represented by the semi-transparent 
area. 


much the inverse trend of the spanning trees ratio around the time 
corresponding to the maximum compression of the system. The 
mass segregation on a longer time scale is slightly less than that 
observed for the same value Q = 0 in Fig. [T] For values of the 
initial virial ratio Q > 0.1 the situation deeply changes. Actually, 
the presence of gas reduces the role of close encounters between 
stars and the process of mass segregation is less efficient. Although 
this point deserves a deeper investigation, it is possible to point out 
that the case Q = 0.1 is profoundly different respect to the case of 
abscence of gas as it is seen also in the evolution of the distribution 
of the velocities of the stars (Fig.|^. 

Figure 1^ shows the distribution of the velocities of the stars 
at different times. In both the cases of presence and of absence of 
a gas component, the initial velocity distributions follows approx- 
imatively the same trend. The situation deeply changes after ~1 
crossing time: the velocity dispersion in the system which does not 
include gas becomes rapidly broader than the other case. This im¬ 
plies that the background potential reduces the efficiency of the 2- 
body interactions and, then, reduces the possibility to form gravity- 
driven substructures (clumps) and the efficiency of mass segrega¬ 
tion. At the bounce, the difference between the two systems is even 
more pronounced and, after ~2 crossing times, the two distribu¬ 
tions are completely different. In particular, the presence of the gas 
tends to pack the velocity distribution toward the same, small, ve¬ 
locity. On the other hand, the velocity distribution, when no gas is 
present, is extended over a huge range of velocities. This is the nat¬ 
ural consequence of more efficient 2-body encounters which make 
the system segregate masses more rapidly. 

We have also studied the case of an expanding background 
gas, represented as a Plummer sphere where the core radius, 
is not constant but, rather, evolving according to the law rdt) = 
rc(0) exp(f/T), where the time scale r is set equal to the crossing 
time of the system {tc) and rc(0) = R. The total amount of the 
mass in gas is given by the star formation efficiency parameter S 
(see Eq.|^; we recall that, in our simulations, we have S — 0.3. 
The law is such that the core radius amplify its value of a factor k 
after a time tk = r In k. Fig.j^gives the spanning trees ratio in this 
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Figure 7. The star velocity distribution in the case of of a system without 
any background gas (Case I, solid line) and in presence of a background 
static gas (dashed line). The initial value of Q is 0.1. The four panels rep¬ 
resent different times: the initial state, the situation after 1 crossing time, at 
the time of the bounce and after 2 crossing times. 
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Figure 8. As in Fig. ^ but with the inclusion of a relic gas treated as a 
smooth gravitational potential in expansion according to a given law. The 
different values of the initial Q are labelled. Each curve represents an av¬ 
erage value of Ajj /taking into account the results obtained from the 
single runs, while the error (standard deviation) is represented by the semi¬ 
transparent area. 

expanding case, for 4 values of the initial virial ratio Q. We note 
how the spanning trees ratio do not show significant differences 
from the behaviour they have in the static gas potential case (Fig. 
1^. In particular, the curves in Fig. [^reproduce the trends reported 
in Fig.|^(the static gas case) for values 0.1 Q ^ 1. Compared to 
the static gas case, here the spanning trees ratio are comparatively 
larger, stabilizing at a value of ~ 1.5 after ~ 8fc, in the zero initial 
velocities case (Q = 0). This means that, also in this case, the 
violent mass segregation effect is quite evident. 


4 DISCUSSION AND CONCLUSIONS 

The main scope of this paper was the study of the topic of quick 
mass segregation in a stellar system composed by a moderate num¬ 
ber of stars. The relevance of this topic is suggested by observations 
showing, unexpectedly, mass segregation in young clusters, often 
still embedded in the gas left by the gaseous protocloud. Despite 
some effort, both on the observational and theoretical side, no firm 
conclusions on the origin of the observed mass segregation at the 
level observed in young embedded clusters have been drawn on a 
theoretical and/or numerical side. 

To attack this problem properly, and without exceeding in 
useless details which can mask the main physical processes gov¬ 
erning the phenomenon, we considered clusters of moderate size 
(128 ^ W ^ 1,024) as composed by stars with two different 
masses (m; and m^, with mu = 2mi), spanning a range of initial 
virial ratios Q between cold conditions, <5 = 0, through cool, and 
warm conditions, Q ^ 1. 

Numerical simulations of the W-body evolutions were per¬ 
formed by mean of the high-precision and high-speed Hi GPUs 
code ( [Capuzzo-Dolcetta et al.|[2013| ). It is worth noting that, al¬ 
though the code, at the moment, does not implement any special 
treatment to integrate neither binaries nor multiples and to deal ac¬ 
curately with very close encounters, the simulations we performed 
so far are numerically very precise despite the criticality of the ini¬ 
tial conditions of the tested systems (initially null velocities of the 
stars, corresponding to the most violent collapse). To quantify this. 


we show in Fig. the relative variation of the total energy of the 
system in the case N — 1, 024, Q — 0, without presence of gas 
but with the presence of the central heavier particle. It is see in the 
Figure that the relative variation of energy, \E{t) — i5(0)|/|i5(0)|, 
is always below 10“® even after the core collapse of the system, 
that occurs at time tec ~ Itc, after which very tight multiple sys¬ 
tems of stars form at the center, leading to a deep increase of the 
total energy of the system. This increase is mainly due to the sig¬ 
nificant reduction of the time step used in the integration of motion 
which yields to a deep increase of the number of iterations per time 
unit, and so a significant increase in the round-off error. 

We gave statistical relevance to our results by making the av¬ 
erage of a set of sampled initial conditions for any of the input 
parameters chosen {N, Q, presence or absence of a black hole and 
of a background gaseous component). Moreover, we quantified the 
level of mass segregation by the use of the powerful minimum span¬ 
ning tree method ( [Allison et al.|200^ coupled to the often used la- 
grangian radii. We have also taken into account the role of escapers 
on our results, accurately treating the escapers as those stars hav¬ 
ing positive individual mechanical energy and far enough from the 
cluster such to have almost zero probability to be recaptured by the 
cluster itself. 

We found that clusters with cool initial conditions (Q < 0.3) 
segregate masses, in a way almost independent of N in the range 
of N studied. This happens in two distinct dynamical phases which 
both contribute to the mass segregation on short or even very short 
time scales. They are: 

(i) the formation of substructures (clumps) as the system un¬ 
dergoes its initial collapse, which segregate masses very quickly 
(this phenomenon acts in the first evolution of the system before 
the bounce of the system while, during the bounce, the clumps are 
completely removed); 

(ii) the subsequent formation of a very dense core which is re¬ 
sponsible for the mass segregation of the system that is clearly evi¬ 
dent after its bounce when no clumpy structures are formed. 

The first phase practically corresponds to the initial gravita¬ 
tional instability of the system and thus is related to the value of 
the Jeans mass of the system in dependence on the assumed initial 
virial ratio, in the sense that warmer initial conditions tend to cancel 
the role of this phase. On the other side, the second phase, although 
it is stronger in its effects for lower values of Q, remains signifi¬ 
cant in segregating masses on an extended interval of Q. Both the 
phases occur on time scales much shorter that the (ill defined, when 
out of equilibrium) 2-body, gentle relaxation time of the system. 

On both the two phases of mass segregation has a relevance the 
inclusion of a “background” potential which mimics the presence 
of a residual gas left after star formation as well as the inclusion 
of an additional body as heavy central object (a stellar black hole). 
Quite clearly, the role on the cluster dynamics played by the back¬ 
ground gas general potential depends on the characteristics of its 
time evolution (if it is very rapid it may influence both the phases). 
The cases studied here (static and expanding gas) indicate that the 
addition of gravitational binding energy tends to smear out the local 
processes that give mass segregation, although segregation is still 
clearly visible. 

The presence of a black hole {rriBH = 25m/i = 50m;) influ¬ 
ences just the second phase and the following, secular, evolution of 
the system. The black hole acts both in reducing the efficiency of 
close encounters between stars, thus decreasing the rate of energy 
exchanges and, so, the resulting mass segregation. 

To conclude: the results presented here are rough if considered 
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Figure 9. The relative energy variation for the most critical case simulated: 
N = 1,024, Q = 0, the inclusion of a central stellar mass black hole and 
without a gaseous component 

on a pure astrophysical side, in what we did not pretend to represent 
in detail real open clusters, but they are interesting enough on the 
side of dynamics of small- to intermediate- A'^-body systems. 

Needed future steps to assess more firmly on the modes of vio¬ 
lent mass segregation would require a proper inclusion of a primor¬ 
dial binary population, of a realistic mass function, a self consis¬ 
tent treatment of the residual gas and of the gas expelled via stellar 
winds during stars’ evolution. Our future aims include also a very 
accurate treatment of close encounters between stars, employing a 
regularization method for the A^-body problem. To do this, we will 
use an updated version of HiGPUs, named HiGPUs-R, which is 
still under development. 
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